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ABSTRACT 

The  U.S.  Navy  has  an  interest  in  the  use  of  laser  systems  for  surface  ships.  Such  systems  must  operate  within 
a  thin  near-surface  environment  called  the  marine  atmospheric  surface  layer.  There  exist  substantial  gradients 
in  temperature  and  momentum  within  this  layer  which  make  turbulence  a  strong  function  of  height.  We  are 
interested  in  robust  and  simple  optical  turbulence  models  that  can  be  used  to  predict  turbulence  along  near¬ 
horizontal  paths.  We  discuss  several  different  models  that  are  based  upon  similarity  theory,  and  we  compare 
the  models  with  field  transmission  data  taken  from  both  over-water  and  over-land  propagation  paths. 

1.  INTRODUCTION 

There  are  several  envisioned  uses  of  laser  systems  within  the  U.S.  Navy.  The  high  bandwidth  capacity  of  point- 
to-point  laser  communications  systems  is  an  appealing  concept.  Similarly  the  use  of  a  high  energy  laser  for 
ship  defense  promises  great  versatility,  since  such  a  system  permits  speed-of-light  delivery  to  a  target.  It  also 
permits  graduated  engagements  that  could  range  from  a  simple  disruption  of  activity  to  destruction. 

In  this  paper  we  discuss  the  analysis  of  the  performance  of  a  laser  system  within  the  marine  atmospheric 
surface  layer.  Determination  of  the  refractive  index  structure  parameter  is  the  most  important  predictor 
of  the  effects  of  optical  turbulence  on  beam  quality,  and  we  will  focus  on  a  comparison  between  models  and 
measurements  of  C%- 

Almost  all  modern  models  for  turbulence  prediction  are  based  upon  the  Monin-Obukhov  similarity  theory. 
We  will  examine  and  critique  two  similarity-based  models  for  C The  primary  method  for  examination  and 
critique  of  the  model  attributes  is  a  direct  comparison  of  two  different  models  against  field  data. 

The  first  data  set  considered  is  taken  from  a  propagation  test  in  the  desert.  This  data  simplifies  the 
examination  of  the  behavior  of  the  model  since  each  diurnal  period  includes  strong  stable  and  strong  unstable 
components,  and  a  very  characteristic  transition  through  neutral  stability.  Furthermore,  aerosol  attenuation  is 
not  a  factor  in  most  of  the  high  desert  data. 

The  second  data  set  considered  is  taken  from  propagation  over  the  ocean  surface.  For  this  data  no  prominent 
diurnal  period  appears  in  the  data  set,  and  both  heat  capacity  and  mixing  are  dominant  factors. 

2.  FIELD  TEST  EQUIPMENT  AND  TEST  PROTOCOL 

The  field  data  examined  in  this  test  was  generated  by  a  transmissometer-scintillometer  comprising  a  blackbody 
source,  a  transmitting  telescope,  a  receiving  telescope  and  both  a  mid-wave  detector  (3.5/rm-4.1/nn)  and  a 
long-wave  (9.8/iin-l  1 .3 /./in)  detector.  The  two  telescopes  are  nearly  identical  with  f/4  20  cm  primary  mirrors 
and  hence  80  cm  focal  lengths.  The  blackbody  is  operated  at  1200°K  with  a  1  KHz  chopper  and  a  radio 
frequency  (162.1  MHz)  reference  signal.  The  blackbody  aperture  is  6.4  mm  yielding  a  field  beam  divergence  of 
ss  8  milliradians. 
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The  receiver  system  consisted  of  the  telescope  and  both  a  long-wave  infrared  detector  and  a  mid-wave 
infrared  detector.  The  signal  from  the  detectors  was  separated  from  the  chopped  carrier  waveform  by  means  of 
a  lock-in  amplifier  system.  The  scintillation  measurement  procedure  sampled  and  recorded  the  signal  at  a  300 
Hz  rate  for  109  seconds.  This  measurement  was  performed  every  15  minutes,  24  hours  a  day. 

3.  THE  CHINA  LAKE  FIELD  TEST:  A  TEST  OVER  LAND 

A  first  set  of  turbulence  data  to  be  examined  here  is  from  a  9  day  test  period  from  12  July  2001  through  20 
July  2001  at  the  Naval  Air  Warfare  Center  at  China  Lake,  CA.  In  spite  of  the  name  of  the  test  location,  the 
test  range  consists  of  a  relatively  homogeneous  flat  and  dry  lake-bed.  These  dates  correspond  to  our  time  scale 
using  day-of-year  193  through  201.  The  transmission  path  was  a  near-horizontal  slant-path,  with  the  receiver 
at  a  height  of  22  m,  and  the  transmitter  height  was  2  m.  The  path  length  for  the  data  used  in  this  paper  was 
3850  meters  in  range. 

There  were  four  distinct  sites  that  served  as  sources  for  meteorological  data  for  the  test  period.  Two  of  the 
sources  provided  data  over  the  full  period  of  the  test,  and  we  have  compared  the  two  data  sets.  One  source 
(‘Tower  8’)  is  approximately  10  km  from  the  propagation  path,  and  the  second  (‘Air  Field’)  is  approximately  3 
km  distant.  This  is  not  an  ideal  configuration  since  one  would  hope  for  co-location  of  the  meteorological  and 
propagation  measurements,  but  the  homogeneity  of  the  area  (a  large  dry  lake  bed)  and  an  inter-comparison  of 
the  data  indicates  that  the  larger  scale  events  recorded  in  both  meteorological  data  records  is  close  to  the  actual 
propagation  conditions.  The  ground  surface  temperature  was  recorded  by  an  infrared  radiometer  located  with 
100  m  of  the  receiver  site. 

The  meteorological  data  is  recorded  at  5  minute  intervals.  We  use  either  3  point  or  5  point  averages  for 
the  meteorological  data  since  the  scintillation  data  is  recorded  at  15  minute  intervals,  and  the  synchronization 
between  the  meteorological  clock  and  the  scintillation  clock  is  expected  to  have  possible  errors  of  several  minutes. 

4.  COMPARISON 

We  chose  data  from  the  China  Lake  site  because  of  the  prominent  meteorological  features  that  appear:  large 
daily  solar  insolation,  and  a  large  diurnal  change  in  wind  speed.  This  creates  a  pattern  in  the  most  important 
model  input  data  and  thus  helps  to  isolate  and  partially  de-couple  some  of  the  complex  features  in  the  model 
equations.  The  over-land  data  shows  the  diurnal  solar  insolation  effects  prominently,  and  in  addition  includes 
strong  diurnal  wind  speed  variations.  The  effects  of  the  wind  and  the  solar  insolation  can  be  distinguished,  and 
both  are  apparent  in  the  C2  calculation. 

Scintillation  data  was  collected  by  our  scintillometer  along  a  shorter  (1280  m)  path,  and  by  a  saturation- 
resistant  commercial  scintillometer  (made  by  Scintec),  which  provided  unsaturated  values  of  C2  such  that 
max{C2}  >  10-13m-2/3.  The  extent  of  the  longer  (3850  m)  propagation  path  that  is  the  source  of  the  data 
investigated  here,  caused  the  C2  data  (plotted  with  squares  in  fig.  1)  to  show  the  effects  of  saturation;  the 
mid-day  values  of  C2  <  10-13m-2/3.  In  the  comparison  of  model  data  with  experimental  data  the  first  point 
to  emphasize  is  the  much  larger  variance  of  the  model  C2  (shown  as  ‘+’  signs)  than  of  the  field  C2  data.  Note 
that  the  large  excursions  of  the  C2  prediction  can  exceed  10-12m-2/3,  and  these  values  are  well  in  excess  of 
any  of  the  measured  maxima  of  C2. 

We  will  use  the  data  for  the  48-hour  period  of  days  195-196  to  examine  the  elements  of  the  model.  There 
is  a  fairly  clear  pattern  to  be  seen  in  the  divergence  of  the  model  and  held  data.  The  first  feature  to  notice  is 
the  occurrence  of  minima  in  C2  occurring  «  30  minutes  after  sunrise,  and  «  30  minutes  before  sunset.  The 
noteworthy  features  in  fig.  1  are:  1)  the  held  data  (black  squares)  show  the  minima  with  the  morning  minima 
less  strong  than  the  evening,  and  the  PAMELA  model  produces  a  stronger  minimum  at  both  points;  2)  the 
large  fluctuations  in  the  modelled  C2  for  the  hrst  half  of  the  day  (from  195.0  to  195.5,  and  from  196.0  to  196.5) 
that  generate  unrealistically  large  C2  values;  3)  a  sustained  and  smooth,  steady  signal  for  the  second  half  of 
each  day  that  is  slightly  under-predicting  the  held  C2  values. 
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Figure  1.  A  48-hour  time  series  from  transmission  across  the  3850  m  path.  determined  by  direct  transmission 
measurements  at  3.85/rm  is  shown  with  the  boxes,  and  PAMELA  model  data  is  shown  with  plus-signs.  Note  the 
prominent  minima  in  C„  occurring  «  30  minutes  after  sunrise,  and  «  30  minutes  before  sunset.  The  upper  panel  shows 
for  wind-speed  threshold  0.1  m/s,  while  the  lower  panel  shows  C„  calculated  for  wind-speed  =  max{uo,2.5}. 

5.  PAMELA  MODEL  DETAILS 

Solar  irradiance  R  is  calculated  using  the  geographical  and  temporal  inputs,  via  spherical  astronomy  calculations. 
The  estimated  solar  irradiance  R  is  used  to  determine  the  radiation  class  cr  =  R/ 300.  For  a  wind-speed  t)o 
define  the  wind-speed  class  cw  =  min{i>o/2,4}  and  then  the  Pasquill  stability  category  P  can  be  determined: 

p  =  ~(4  ~cw+cr) 

2  1  } 

The  roughness  length  for  the  shrub-covered  dry  lakebed  can  be  estimated  from  tables  such  as  that  found  in 
Arya1:  zr  =  0.16m,  and  from  this  it  is  possible  to  calculate  the  Obukhov  buoyancy  length  scale  L : 
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where 


L  = 


(ai  P  +  a2P3)^-(a3"04|p|+05p2) 


(2) 


ai  =  0.004349,  a2  =  0.003724, 
a3  =  0.5034,  a4  =  0.231, 
a5  =  0.0325. 

The  eddy  viscosity  (diffusivity)  relations  are  defined  next.  For  mean  vertical  velocity  W  and  fluctuating 
part  w,  mean  horizontal  velocity  U  and  fluctuating  part  u,  define  vertical  momentum  flux  in  terms  of  the  eddy 
viscosity  K.m : 

—  K  -'dU " 
uw  =  — 


dz 


(3) 


For  the  mean  potential  temperature  0  and  fluctuating  part  0,  define  vertical  heat  flux  in  terms  of  the  eddy 
diffusivity  of  heat  Kft: 

0^  =  -Kh  ( ^ 


dz 


(4) 


For  the  mean  specific  humidity  Q  and  fluctuating  part  q  define  the  vertical  water  vapor  flux  using  the  eddy 
diffusivity  of  water  vapor  Kw : 

—  IS  (dQ" 
qw  =  —Kw  — 


dz 


(5) 


The  dimensionless  wind  shear  ',)  and  the  dimensionless  potential  temperature  gradient  <ph{C)  are  ex¬ 
pressed  as  functions  of  the  scaled  buoyancy  parameter  £  =  z/L.  The  turbulent  exchange  coefficients  for  heat 
Kh  and  momentum  Km  are: 


KU*Z  KU*Z 

=  7TT  cXnCl  i\-rn,  = 


4>h{  0 


(6) 


where  k  —  0.4  is  the  von  Karman  constant.  Note  that  in2  Ki,  =  Km.  From  the  wind  speed  vo  and  the  roughness 
length  zr  the  friction  velocity  it*  is  given  by: 


u*  = 


KV0 


and  the  characteristic  temperature  T* 


T*  = 


ln(z/zr) 
-H 


Cppu* 


(7) 


(8) 


is  defined  using  the  heat  flux  H,  specific  heat  cp,  and  mass  density  p.  With  pressure  Pa  measured  in  millibars, 
the  atmospheric  refractive  index  n: 


n  —  1 

The  field  data  uses  A  =  3.85/jm,  and 

The  eddy  dissipation  rate  e  is: 


77.6  x  10-6Po 
T 


1  + 


7.52  x  10 
A2" 


-3 


(9) 


dn  77.6  x  10  6PaT*(/)h 
~dz  ~  OAzT2 


e  = 


^*(^m  0 

0.42 


(10) 
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Figure  2.  A  comparison  of  T*,  tt*,  and  for  a  24-hour  period.  For  this  data,  the  minimum  wind  speed  threshold  is 
set  to  0.1  m/s. 


and  with  the  constant  b  ss  2.8  we  can  evaluate  C%: 


2  _  bKh  / dn\2 
n  ~  e1/3  \dz) 


(11) 


To  elucidate  these  features,  in  fig.  2  the  important  elements  T*,  it*,  and  C2n  are  compared.  To  facilitate  the 
display  of  the  graphs  of  the  three  quantities,  the  time  domain  for  the  comparison  has  been  restricted  to  a 
24-hour  period:  195.0-196.0. 

As  noted  in  the  figure  legend,  the  curve  is  shown  as  log(C^),  and  it  has  been  shifted  to  overlay  the  graphs 
of  T*  and  it».  As  can  be  seen  from  eqns.  (7)  and  (13),  — »■  oo  as  wind-speed  Vo  — >  0,  so  the  minimum  wind 

speed  must  be  bounded  away  from  zero.  In  fig.  2  the  wind-speed  threshold  is  0.1  m/s.  In  this  plot,  we  averaged 
the  C„  calculation  over  a  25-minute  span  to  damp  some  of  the  more  extreme  spikes.  This  averaging  was  done 
simply  to  reduce  the  sharp  fluctuations  in  the  T*  and  C%  graphs.  It  is  still  the  case  that  the  C%  calculation  is 
highly  fluctuating. 

The  other  two  curves  are  included  to  elucidate  the  factors  in  the  behavior  of  the  calculation.  First,  note 
that  eqn.  (7)  implies  that  the  fluctuations  in  the  graph  of  u*  are  essentially  showing  the  variations  of  wind  speed 
vo-  As  the  wind  speed  oscillates  around  a  small  value  (the  time  from  195.0  to  195.3),  there  are  large  fluctuations 
induced  in  both  T*  and  in  C^-  The  direct  relationship  can  be  seen  by  expanding  eqn.  (11)  as: 


rj-iZ 


=  1.84  b(j>h 


1  \  1/3  /  77.6  x  10-6Pa 


(</>m  -  <) 


J'  2 


,-2/3 


-H 

Cppu, 


(12) 

(13) 


For  this  model  experiment,  we  hold  the  height  z  =  10m.  The  important  variations  in  the  signals  and  T* 
are  thus  entirely  produced  by  the  fluctuations  of  the  values  of  heat  flux  H  and  u*  which  (see  eqn.  (7))  is  just  a 
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linearly  scaled  wind  speed  vq.  In  the  examination  of  the  model  behavior  in  the  following  sections,  we  will  follow 
the  values  of  the  variables  it*,  T*,  and  H  which  are  indicated  by  the  under-braces  in  eqns.  (12)  and  (13). 

The  point  at  which  heat  flux  H  changes  sign  is  evident  in  fig.  2,  since  T*  changes  sign  at  the  same  point. 
The  effect  of  small  values  of  it*  is  also  apparent  for  195  <  time  <  195.5:  as  u*  becomes  small,  both  T*  and  C\ 
experience  large  spikes  in  value.  Around  noon,  (195.5)  the  mean  wind  speed  starts  to  increase.  The  effect  on 
both  C%  and  T*  is  apparent  in  fig.  3  as  the  large  fluctuations  in  the  graphs  of  both  are  greatly  reduced  for  time 
>  195.5. 


6.  DISCUSSION  OF  THE  PAMELA  MODEL 

The  difficulty  with  the  PAMELA  model  is  that  removal  of  the  sensitivity  to  wind  speed,  when  the  wind  speed 
is  low,  is  not  a  simple  modification.  If  we  insert  a  higher  threshold  for  wind  speed,  as  is  done  in  fig.  3  the  large 
scale  fluctuations  are  still  unrealistically  strong.  In  the  upper  panel  of  fig.  3  the  same  data  is  shown  for  a  wind 
speed  threshold  of  1.0  m/s,  thus  wind  speed  i>o  =  min{t>o,  1.0}  m/s.  In  this  figure,  the  signal  was  also  averaged 
to  smooth  out  the  data. 

The  effects  of  low  wind-speed  on  become  even  more  apparent  when  the  wind  speed  threshold  is  increased 
to  2.5  m/s  in  the  lower  panel  of  fig.  3.  Although  this  threshold  ensures  that  the  extreme  values  of  the  model  C/ 
are  suppressed,  many  of  the  other  legitimate  dynamic  elements  of  the  signal  are  also  removed  by  this  process. 

We  recognize  that  there  are  other  dependencies  that  affect  the  amplitude  and  shape  of  the  prediction. 
We  have  not  considered  for  example  the  contributions  from  dimensionless  wind  shear  <j)m  and  the  potential 
temperature  gradient  (f>h ■  But  these  are  functions  of  the  Monin-Obukhov  length,  and  should  not  show  the 
instantaneous  effects  of  wind  speed. 

7.  BEHAVIOR  OF  PAMELA  FOR  PROPAGATION  OVER  WATER 

In  this  section  we  examine  the  application  of  the  PAMELA  model  to  the  prediction  of  optical  turbulence  over 
the  ocean.  A  data  set  generated  from  measurements  on  a  7  km  path  across  Zuniga  Shoals  outside  of  San 
Diego  Bay  was  used  as  the  data  for  this  test.  A  record  of  scintillation  data  over  a  two-day  period  is  shown  in 
fig.  4.  The  field  test  equipment  for  the  scintillation  measurements  was  quite  similar  to  the  system  described  in 
section  2  above.  A  scintillation  measurement  was  made  every  20  minutes,  24  hours  per  day  throughout  the  test. 
Simultaneous  meteorological  measurements  were  made  from  a  flux  buoy  anchored  near  the  mid-point  of  the 
propagation  path,  including  in  particular  air  temperature,  sea  surface  temperature,  humidity  and  wind-speed 
measurements. 

Note  that  the  field  data  shows  no  evidence  of  the  diurnal  change  in  sign  of  the  surface  heat  flux.  The 
ocean  has  a  very  large  heat  capacity,  and  effective  mixing  processes.  Thus  solar  insolation  is  only  one  of  several 
important  factors  in  the  determination  of  sea  surface  temperature  and  subsequently  Ta;r  —  Tsurface.  In  fact 
Tail-  —  Tsurface  <  0  can  persist  for  days  or  weeks  in  an  open  ocean  environment.  For  the  comparison  shown  in 
fig.  4  we  set  the  PAMELA  parameters  to  a  theoretically  more  benign  range:  cloud  cover  throughout  the  five 
day  period  was  set  to  overcast,  the  roughness  length  z0  =  0.01,  and  we  use  a  wind-speed  threshold  Vq  =  2.5 
m/s. 

The  effect  of  full  cloud-cover  in  the  PAMELA  model  is  to  make  the  mid-day  local  maximum  in  still  less 
than  the  overall  maximum  C^-  However  it  is  apparent  from  fig.  4  that  the  predicted  over  the  full  48-hour 
period  exceeds  the  measured  values  by  one  to  two  orders  of  magnitude.  Clearly  the  dynamic  range  of  PAMELA 
simulation  is  excessive  for  data  that  is  generated  in  over-water  environments,  especially  considering  that  a  full 
cloud-cover  was  invoked. 

The  PAMELA  model  predicts  that  the  zero-flux  conditions  near  sunrise  and  sunset  will  occur,  but  this  is 
not  seen  in  the  data.  This  points  out  a  weakness  in  the  current  PAMELA  model:  characteristics  such  as  heat 
flux  at  the  air-sea  interface  cannot  be  represented.  The  dynamical  features  in  the  graph  of  the  PAMELA  model 
in  fig.  4  are  almost  entirely  missing  in  the  graph  of  the  recorded  infrared  data,  and  in  fact  a  cross-correlation 
calculation  between  the  model  log(C^)  and  the  measured  transmission  log(C^)  confirms  that  the  two  time  series 
are  essentially  uncorrelated. 
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Figure  3.  A  comparison  of  T*,  «*,  and  C2.  For  the  upper  panel,  the  minimum  wind  speed  threshold  is  set  to  1.0  m/s. 
In  the  lower  panel,  the  minimum  wind  speed  threshold  is  set  to  2.5  m/s. 

8.  THE  NAVY  SURFACE  LAYER  OPTICAL  TURBULENCE  MODEL-NSLOT 

The  refractive  index  structure  parameter  C2  is  defined  implicitly  by  the  refractive  index  structure  function: 

[■ n{x )  —  n{x  +  r)]2  =  C2r2^3  (14) 

where  the  turbulence  is  assumed  to  be  isotropic  and  hence  dependent  only  on  r  =  |r|. 

Optical  turbulence  is  generated  by  fluctuations  in  the  refractive  index  n.  The  derivation  of  Andreas3  defines 
refractive  index  fluctuations  in  terms  of  fluctuations  in  temperature  T  and  humidity  q. 

n'  =  A( A,  P,  T,  q)T'  +  B( A,  P,  P,  q)q'  (15) 


Threshhold  wind-1.0  m/s 
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Figure  4.  A  comparison  of  C2  over  the  coastal  Pacific  Ocean  for  days  253  through  255.  Overcast  skies  (complete  cloud 
cover),  a  roughness  length  zT  =  0.01  and  a  threshold  wind-speed  Vo  =  2.5  m/s  were  used  for  the  PAMELA  model. 


where  a  prime  indicates  an  instantaneous  turbulent  fluctuation  of  a  quantity  about  its  mean  value,  and 


A  =  ^  =  -10  jrni(A)  +  [m2( A)  -  mi  (A)]  ^  j 

B=-  =  1°-6[m2(A)-mi(A)]^ 

where  7=  1  +  0.6078(7,  e  is  the  ratio  of  ideal  gas  constants  for  dry  air  over  water,  and 


mr  (A)  =  23.7134  + 


6839.397 
130  -  A-2 


45.473 
38.9  —  A~2 


to2(A)  =  64.8731  +  0.58058A-2  -  0.007115A"4  +  0.0008851A"6 


(16) 

(17) 


(18) 

(19) 


with  A  given  in  /rm. 

With  the  structure  parameter  for  temperature  C2 ,  the  humidity  structure  parameter  C2,  and  the  temperature- 
humidity  structure  parameter  Ctq,  we  follow  Andreas3  and  define  : 


C2  =  A2  Cl  +  2  ABCtq  +  B2C2q 


A  method  to  estimate  C2  from  routine  meteorological  observations  is  available  in  the  Navy  Surface  Layer  Optical 
Turbulence  (NSLOT).4  Let  the  A-operator  denote  a  mean  air-sea  difference  for  a  particular  meteorological 
measurement.  Thus  an  ‘air’  measurement  can  be  made  within  the  surface  layer,  a  small  distance  above  the  sea 
surface,  and  a  second  ‘sea’  measurement  taken  at  the  surface  to  provide  the  two  measurements  indicated  by  the 
A-operator. 


,2  f(£)k2(A2AT2  +2ABrr^ATAq+B2Aq2) 

n~  z^Mz/z0T)-^Tm 


(20) 
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Figure  5.  The  upper  panel  shows  a  comparison  of  C„  calculation  with  mid-wave  infrared  propagation  data  over  Zuniga 
Shoals  (a  7  km  over-water  path)  shown  with  squares,  and  the  NSLOT  model  prediction  shown  with  signs.  The  lower 
panel  follows  AT  =  Ta ir  —  Tsurface  for  the  same  time  period. 


where  zqt  is  the  temperature  roughness  length,  rxq  is  the  temperature-specific  humidity  correlation  coefficient, 
\E't(0  is  a  dimensionless  profile  function,  and 


Hi)  =  MO  =  M(0  =  /,({)  =  {  +  f ;  °  (21) 

where 

,  =  zg(AT  +  Q.QlTAg)[hi{z/z0U)  -  ^r/(g)]2  , ^ 

4  evAU2Mz/z0T)-^Tm 

The  behavior  of  the  NSLOT  model  for  data  generated  over  water  is  shown  in  fig.  5.  In  the  figure,  the  values 
of  field  observations  of  are  compared  to  predictions  of  the  NSLOT  model  parameterized  from  meteorological 
data  taken  from  a  flux  meteorological  buoy  anchored  at  the  mid-point  of  the  propagation  path.  The  lower 
panel  in  the  figure  shows  the  values  of  AT  =  Tajr  —  Tsurface  over  the  analyzed  period,  and  for  this  period, 
AT  =  Ta jr  —  Tsurface  >  0,  and  the  NSLOT  model  over-predicts  CZ  values.  It  can  be  seen  that  predictions  by 
NSLOT  (shown  as  small  +  signs)  are  in  general  larger  C\  than  the  measured  values  (seen  as  plotted  squares). 
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Figure  6.  The  upper  panel  shows  a  comparison  of  C„  calculation  with  mid-wave  infrared  propagation  data  over  Zuniga 
Shoals  shown  with  squares,  and  the  NSLOT  model  prediction  is  shown  with  “+’  signs.  The  lower  panel  follows  AT  = 
Tail-  —  Taurface  for  the  same  time  period.  Data  points  satisfying  —0.1  <  AT  =  Tair  —  Tsurface  <  0.2  are  over-plotted  with 
a  diamond  in  both  panels. 

In  figure  6  is  shown  a  comparison  of  NSLOT  predictions  with  scintillation  measurements  at  a  different  period 
of  2.5  days  earlier  in  the  same  test.  For  the  first  day  (253.0-254.0),  AT  =  Tajr  —  Tsurface  «  0,  and  this  causes 
NSLOT  to  predict  excessively  small  values  of  C 2.  This  behavior  can  be  seen  in  eqn.(20)  by  noting  that  the  two 
terms  highlighted  with  over-brackets  go  to  zero  as  AT  —>  0. 

To  illuminate  in  fig.  6  how  the  near-zero  values  generate  poor  C2  estimates  in  NSLOT,  the  points  satisfying 
—0.1  <  AT  =  Ta;r  —  Tsurface  <  0.2  are  over-plotted  with  a  diamond  both  in  the  lower  panel  and  in  the  upper 
panel.  Those  points  near  zero  correspond  to  10-17m-2/3  <  C2  <  10-16m-2/3,  and  these  points  deviate  from 
the  measured  values  in  some  cases  by  more  than  a  factor  of  ten.  The  general  trend  is  for  the  NSLOT  model  to 
under-predict  the  measured  C2  values  for  AT  near  zero. 

In  the  remaining  1.5  days,  (254.0-255.5)  the  NSLOT  model  does  quite  well  at  predicting  the  observed  C2 
values.  Thus  NSLOT  is  successful  at  prediction  when  the  conditions  are  in  a  condition  common  to  the  open- 
ocean:  AT  =  Tajr— Tsurface  <  0.  This  provides  a  striking  contrast  with  the  PAMELA  model  which  was  compared 
to  the  same  data  (253.0-255.0)  as  shown  in  fig.  4. 
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9.  DISCUSSION  OF  THE  NSLOT  MODEL 


As  noted  above  the  NSLOT  model  performs  well  for  a  common  open-ocean  condition:  AT  =  Tajr  —  Tsurface  <  Su 
where  for  example  8U  =  — 0.4°C.  The  model  also  appears  to  over-predict  for  AT  =  Ta;r  —  Tsurface  >  Ss,  where 
for  example  Ss  =  0.6°C.  (The  selections  for  Su  and  5S  are  simple  rough  estimate  guided  by  evidence  from 
Frederickson4).  For  the  in-between  region  Sa  <  AT  =  Tair  —  Tsurface  <  Su  NSLOT  tends  to  under-predict, 
sometimes  by  several  orders  of  magnitude. 


10.  CONCLUSIONS 

We  have  tested  two  surface  layer  optical  turbulence  models,  PAMELA  and  NSLOT,  both  of  which  have  been 
developed  from  surface  layer  similarity  theory.  The  first  model,  PAMELA,  is  a  relatively  simple  model  to 
implement.  However,  there  are  two  problems  that  appeared  in  the  PAMELA  test.  First,  the  model  is  highly 
sensitive  to  low  wind-speeds.  The  second  problem  is  the  lack  of  a  surface  characterization  in  PAMELA,  which 
prevents  a  successful  prediction  for  maritime  conditions.  We  are  unable  to  adjust  parameters  in  the  model  to 
yield  an  accurate  scintillation  prediction  for  marine  propagation  conditions. 

NSLOT  has  been  developed  specifically  for  the  marine  environment,  and  because  of  this  certain  parameter 
values  such  as  surface  roughness  have  been  ‘hard-wired’  into  the  model.  For  this  reason  we  did  not  test  the 
model  against  over-land  data  in  this  report.  In  tests  against  over- water  data  we  found  that  NSLOT  can  badly 
under-predict  for  AT  =  Ta,r  —  Tsurface  «  0.  As  we  indicated  in  the  section  above,  the  source  of  the  difficulty 
can  be  seen  in  over-bracketed  terms  of  eqn.  (20).  As  AT  — >  0,  the  two  terms  containing  AT  become  negligible 
compared  to  B2Aq2.  A  more  detailed  study  of  this  weakness  in  NSLOT  is  presented  in  the  paper  by  Frederickson 
et  al.4 

The  NSLOT  model  is  in  general  a  more  accurate  model  for  marine  propagation  conditions.  In  comparing 
NSLOT  to  PAMELA,  it  is  important  to  note  that  NSLOT  manages  the  small  wind-speed  regime  more  success¬ 
fully  than  PAMELA.  It  is  apparent  that  the  NSLOT  approach  to  capturing  the  surface  heat  flux  by  means  of 
Tair  and  Tsurface  measurements  and  thence  a  direct  determination  of  Tajr  —  Tsurface  is  preferable  to  PAMELA’S 
use  of  a  computed  heat  flux  via  solar  insolation. 

Both  models  are  ultimately  constrained  by  the  limitations  of  surface  layer  similarity  theory.  In  particular, 
the  theory  requires  wind  speeds  that  are  not  too  small,  and  a  surface  layer  of  sufficient  thickness  to  encompass 
the  air  temperature,  humidity,  and  wind-speed  measurements.  Conditions  that  include  a  temperature  inversion 
or  Tair  —  Tsurface  >  0  will  be  problematic  for  both  models. 
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